Control and mixing for 2D Navier-Stokes 
equations with space-time localised noise 

Armen Shirikyan* 
October 5, 2011 

Abstract 

We consider randomly forced 2D Navier-Stokes equations in a bounded 
domain with smooth boundary. It is assumed that the random perturba- 
tion is non-degenerate, and its law is periodic in time and has a support 
localised with respect to space and time. Concerning the unperturbed 
problem, we assume that it is approximately controllable in infinite time 
by an external force whose support is included in that of the random force. 
Under these hypotheses, we prove that the Markov process generated by 
the restriction of solutions to the instants of time proportional to the 
period possesses a unique stationary distribution, which is exponentially 
mixing. The proof is based on a coupling argument, a local controllability 
property of the Navier-Stokes system, an estimate for the total variation 
distance between a measure and its image under a smooth mapping, and 
some classical results from the theory of optimal transport. 
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1 Introduction 

The main results of this paper can be summarised as follows: first, suitable con- 
trollability properties of a non-linear PDE imply the uniqueness and exponential 
mixing for the associated stochastic dynamics and, second, these properties are 
satisfied for 2D Navier-Stokes equations with space-time localised noise. To be 
precise, let us consider from the very beginning the 2D Navier-Stokes system in 
a bounded domain Z? C with smooth boundary dD: 

u+ {u,\7}u~ lyAu + Vp ^ f{t,x), divu = 0, x e D, (1.1) 



Here u — {ui,U2) and p are unknown velocity and pressure of the fluid, v > 
is the viscosity, and / is an external force. Let us assume that / is represented 
as the sum of two functions h and rj, the first of which is a given function that 
is smooth in space and time and has a locally bounded norm, while the 
second is either a control or a random force: 



f{t,x)^h{t,x)+Tj{t,x). (1.3) 

In both cases, we assume that rj is sufficiently smooth and bounded, and its 
restriction to any cylinder of the form Jk y< D with Jk = [fc — 1, k] is localised in 
both space and time (see below for a more precise description of this hypothesis) . 
Let us denote by n the outward unit normal to the boundary dD and introduce 
the space 

H ^{ue L^{D,R^) : divu = in D, {u,n} = on dD}, (1.4) 

which will be endowed with the usual norm || • ||. It is well known that for 
any uq <E H problem (jl.ip . (|1.2p supplemented with the initial condition 

m(0, x) — uq(x) (1-5) 
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has a unique solution u — u(t]UQ, f), which is a continuous function of time 
valued in H. 

To state our first result, we introduce a localisation hypothesis on the con- 
trol T]. Let Q C Ji X be an open set and let Qk = {{t, x) : {t~ k + l,x) Q}. 

Localisation. For any integer fc > 1, the restriction of rj to the cylinder Jk x D 
is supported by Qk- 

The theorem below on the stabilisation to a non-stationary solution for 
Navier-Stokes equations is an immediate consequence of a contraction prop- 
erty established in Theorem 13. II 

Theorem A. For any p > and a > there is a finite- dimensional subspace 
E C H^{Q,^^) and positive constants C and d such that the following asser- 
tions hold for any functions E H and h G -ffjQ^(M-|- x ZJjM^) satisfying the 
inequalities 

\\uo\\<p, \\h\\m(Jt,xD) < P for all k> I. 

(i) For any uq E H satisfying the condition ||uo — uqII ^ d there is a control rj 
such that the restriction of rj{t + k — \,x) to Ji x D belongs to E for any 
integer k>l and 

\\u{t\UQ,h + ri)^u{t]Uo,h)\\<Ce-°'^\\uo~uo\\, t>0. (1.6) 

(ii) The mapping (uq, uq, h) ^ rj is Lipschitz continuous for appropriate norms. 
Moreover, 

h||L^(Q„R2) <Ce-""||«o-iio||, k>l. (1.7) 

Let us note that a similar result was proved earlier in |BRS11| for the 3D 
Navier-Stokes system (see also |Fur04i IBar03i IBT04| for some results on stabil- 
isation to a stationary solution). Here we establish some additional properties 
for the control. Namely, we show that it is finite-dimensional in both space and 
time and can be chosen to be a smooth function with respect to the force and 
the initial state corresponding to the reference solution. 

We now turn to the stochastic problem. Informally speaking, we assume 
that the deterministic force h is 1-periodic in time, and the stochastic force rj 
satisfies the following three conditions. Let us denote by rj^ the restriction of 
rj{t -\- k — l,x) to the domain Ji x D. 

Independence. The functions rjk form a sequence of i.i.d. random variables 
in iJi(Ji X D,R^). 

Non-degeneracy. There is a compact subset /C C Hq{Q,'S?) such that the sup- 
port of the law for rjk coincides with /C, and the law of the projection of rjk 
to any finite-dimensional subspace of L'^{Q, M^) has a regular density with 
respect to the Lebesgue measure. 
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Approximate controllability. There is u G H such that problem (|l.ip - (|1.3p is 
approximately controllable to ii with a control function fj such that for 
any k > 1 the restriction of fj{t + k — 1, x) to Ji x D belongs to K., and the 
time of control can be chosen the same for the initial functions uq from a 
given bounded subset. 

We refer the reader to Section 12.11 for the exact description of the hypotheses 
imposed on the stochastic system p.ip - (ll.3p . Note that the hypothesis on 
independence implies that the restrictions of solutions to the integer times form a 
Markov chain in H. The following theorem describes the long-time asymptotics 
of this chain. 

Theorem B. Under the above hypotheses, the Markov chain associated with 
problem pTTIl -lfO l) has a unique stationary measure fi. Moreover, there are 
positive constants C and 7 such that, for any 1-Lipschitz function F : H W 
and any uq G H , we have 



Let us note that the problem of ergodicity was studied intensively in the 
last fifteen years. First results in this direction were established in |FM95l 
[E5001 IEMS01[rBKL02] . and we refer the reader to the book [E5TT] for further 
references. Most of the results established so far concern the situation in which 
the random force is non-degenerate in all determining modes of the problem. In 
the case when the equation is studied on the torus and the deterministic force 
is zero, it was proved in [HM06, HMOS^ that the Navier-Stokes dynamics is 
exponentially mixing for any v > 0, provided that the noise is white in time and 
has a few non-zero Fourier modes as a function of x (thus it is finite-dimensional 
in space and infinite-dimensional in time). Theorem B applies to all v > 0, and 
the random noise does not act directly on the deterministic modes. However, 
in contrast to [HM06[ IHM08) , our noise is localised not in the Fourier modes 
with respect to x, but in the physical space and in time. Furthermore, our 
approach does not use the Malliavin calculus and is based on a detailed study of 
controllability properties for the Navier-Stokes system and a result on the image 
of measures on a Hilbert space under finite-dimensional transformations. To the 
best of my knowledge. Theorem B provides a first result on mixing properties 
for Navier-Stokes equations with a space-time localised noise. 

In conclusion, let us mention that Theorems A and B remain valid in the 
case when the control and the noise act through the boundary of the domain. 
This situation will be addressed in a subsequent publication. 

The paper is organised as follows. In Section[51 we formulate the main result 
of this paper on exponential mixing for the Navier-Stokes system with space- 
time localised noise and outline its proof. Section [3] is devoted to studying a 
control problem associated with the stochastic system in question. The details 
of proof of the main result are given in Section |4l The appendix gathers some 
auxiliary results used in the main text. 




fc>0. (1.8) 
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Notation 

For an open set Q of a Euclidean space, a closed interval J C K, and Banach 
spaces X CY, we introduce the following function spaces. 

LP = LP{Q) is the Lebesgue space of measurable scalar or vector functions 
on Q whose p**^ power is integrable. We shall sometimes write U'{Q,M.'^) to 
emphasise the range of functions. In the case p = 2, the corresponding norm 
will be denoted by || • ||. 

= H^{Q) is the Sobolev space of order s with the usual norm || • H^. As 
in the previous case, we use the same notation for spaces of scalar and vector 
functions. 

Hq = Hq{Q) is the closure in of the space of infinitely smooth functions 
with compact support. 

Bx{R) stands for the ball in X of radius R centred at zero. 

LP{J,X) is the space of Borel-measurable functions u : J ^ X such that 

i/p 



\lp(j,x) = (^Jju{t)fxd?j <oo; 



in the case p = oo, this norm should be replaced by ||u||oo = ess supj 

W^'P{J, X) is the space of functions u € Lp{J,X) whose derivative belongs to 

LP{J,X). It is endowed with a natural norm. 

W{J, X, Y) is the space of functions u G L^{J, X) such that dtu € J, Y). 

C{X, Y) is the space of continuous linear operators from X to F with the natural 
norm. In the case X = Y, we write C{X). 

Cb{X) stands for the space of bounded continuous functions F : X — M; it is 
endowed with the norm 

= sup \F{u)\. 

Li,{X) denotes the space of functions F e Cb{X) such that 
\\Fh:= ||F||oo+supM!i)^<oo. 

u^v \\U-V\\x 

'P(X) is the set of probability Borel measures on X. The space V{X) is endowed 
with the topology of weak convergence, which is generated by the dual-Lipschitz 
metric 

IImi -/"2||l := sup \{F, fj,i) - {F, ij,2)\, 

where {F, p) stands for the integral of F over X with respect to ji. 

We denote by a bounded domain with boundary dD. For T > 0, we set 
Jt = [0,T] and Dt = Jt x D. The following functional spaces arise in the 
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theory of Navier-Stokes equations: 



V^Hr^Hl{D), Xt^W{Jt,V,V*), 
ys^T = {ueXT. u|(5,T) e W{{S, T), V n H\ V)], 

where 5 G (0,T), H is defined by (|1.4I) . and T^* denotes the dual space of V 
(identified with a quotient space in i?~^(Z?,R^) with the help of the scalar 
product in L^). These spaces are endowed with natural norms. 

2 Main result and scheme of its proof 
2.1 Exponential mixing 

Let _D C be a bounded domain with a C^-smooth boundary dD and let 
Di — Ji X D. Consider the Navier-Stotes system with the Dirichlet 

boundary condition (II. 2p and an external force of the form (jl.3|) . We assume 
that h € i?jQ^(R+ X D,M.'^) is a given function which is 1-periodic in time and rj 
is a stochastic process of the form 

oo 

r]{t,x) = "^Ik{t)f]k{t-k+l,x), t>0, (2.1) 
fe=i 

where Ik is the indicator function of the interval (fc — 1, fc) and {rjk} is a sequence 
of i.i.d. random variables in L^(L)i,M^) that are continued by zero for t ^ Ji. 
It is well known that the Cauchy problem for (|1.1I) - (|1.3I) is globally well posed. 
Namely, for any initial function uq £ H there is a unique random process {u,p) 
whose almost every trajectory satisfies the inclusions 

u e C(M+, ff) n Ll^iR+,V), p e 

Equations (|l.ip . and the initial condition 

u{0,x) ^ua{x), X e D. (2.2) 

In what follows, we shall drop the p component of solutions and write sim- 
ply u{t). Let us denote hy S : H x L'^{Di,R'^) — > H the operator that takes a 
pair of functions (mq, /) to u{l), where u{t) is the solution of (jl.ll) . (|1.2p . (|2.2p . 
Well-known properties of 2D Navier-Stokes equations imply that is a contin- 
uous mapping. Moreover, the range of S is contained in V, and the mapping 
S : H X L^(Di,R^) — > F is uniformly Lipschitz continuous on bounded subsets; 
e.g., see Chapter III in |Tem79l . 

Let us consider a solution u{t) of (|l.ip ~ (|1.3p and denote ut — u{k). What 
has been said implies that 

Uk ^ S{uk-i,h + i]k), fc > 1. (2.3) 

Since rjk are i.i.d. random variables in i^(L)i,M^), Eq. (|2.3p defines a homoge- 
neous family of Markov chains in H, which is denoted by {uk,Fu), u £ H. Let 
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Pk{u,r) be the transition function for the family {uk,'Pu) and let : Cb{H) — >■ 
Cb{H) and : V{H) 'P{H) be the corresponding Markov semigroups. 

Let us fix an open set Q C Di and denote by {(fj} C H^{Q,R'^) an or- 
thonormal basis in L'^{Q,R'^). Let x G C'^(Q) be a non-zero function and let 
ipj = Xf]- In what follows, we shall assume that {V'j} are linearly independenl(3, 
and the function h and random process rj satisfy hypotheses (HI) and (H2) for- 
mulated below. Note that if {^Pj} is a complete set of eigenfunctions of the 
Dirichlet Laplacian in Q, then the functions V'j are linearly independent for 
any choice of x- This is an immediate consequence of the unique continuation 
property of solutions for elliptic equations; see Theorem 8.9.1 in | Hor63 |. 

(HI) Structure of the noise. The random variables rjk can be represented 
in the form 

oo 

77fe(t, x)=Y^ bj^jki^jit, x), (2.4) 
i=i 

where ^jk are independent scalar random variables such that \^jk\ < 1 
with probability 1, and {bj} C M is a non- negative sequence such that 

oo 

B:^Y.^,yj,h<^. (2.5) 

Moreover, the law of ^jk possesses a C^-smooth density pj with respect to 
the Lebesgue measure on the real line. 

Let us denote by /C C L^{Q, K^) the support of the law of rjk- The hypotheses 
imposed on rjk imply that is a compact subset in Hq{Q,R'^). Continuing 
the elements of IC by zero outside Q, we may regard as a compact subset 
of i/l(^l,K')• 
(H2) Approximate controllability. There is u G iJ such that for any pos- 
itive constants R and e one can find an integer / > 1 with the following 
property: given v G Bh{R), one can find (^i, . . . ,(^i G /C such that 

||5,(«,Ci,...,0)-^ll <£, (2.6) 

where Si{vXi, ■ ■ ■ ,Ci) stands for the vector ui defined by (12.31) with r/k = Cfc 
and uq — V. 

The following theorem, which is the main result of this paper, establishes the 
uniqueness and exponential mixing of a stationary distribution for the Markov 
family generated by (|2.3p . 

^The assumption on the linear independence of {tpj} is not really needed, and Theorem 1 2. II 
below remains true without it. We make, however, this assumption to simplify the proof of 
Proposition 12.51 



7 



Theorem 2.1. Assume that h G x is 1-periodic in time, and 

Conditions (HI) and (H2) are satisfied. In this case, there is an integer N > I, 
depending on ||^||_f/i(Ui), B, and v, such that if 

b,^0 for j^l,...,N, (2.7) 

then the following assertions hold. 

Existence and uniqueness. The Markov family (ufc,P,i) has a unique sta- 
tionary distribution /i G V{H). 

Exponential mixing. There are positive constants C and 7 such that 

\\Pkiu,-)~ ^i\\l<C{l + \\u\\)e~^'' for allueH, k>0. (2.8) 

Note that condition (|2.7p expresses the space-time non-degeneracy of the 
noise. Thus, the property of exponential mixing holds true even for noises 
whose space-time dimension is finite. 

The general scheme of the proof of Theorem 12.11 is outlined in Section 12.21 
and the details are given in Sections [3] and |4l Here we consider two examples 
for which Condition (H2) is fulfilled. 

Example 2.2. We claim that there is (5 > such that if the 1-periodic function h 
satisfies the inequality ||/i||i2(£i^) < S, then Condition (H2) is fulfilled, provided 
that /C contains the zero element. 

Indeed, if ||/i||i2(2?i) is sufficiently small, then problem (|l.ip . (I1.2p with f — h 
has a unique solution u{t, x) defined throughout the real line and 1-periodic in 
time. To see this, it suffices to take a sequence {un} of solutions for the problem 
in question such that m„(— rt) = and to prove that it converges as rt ^ 00 in the 
space W{[—N, N], V, V*) for any N > 0. The limiting function u is the required 
1-periodic solution. Using standard estimates for Navier-Stokes equations, it is 
easy to prove that 

\\u\\wiJi,v,V') < c{S) ^ as (5-^0. 

This implies that u is globally exponentially stable as t — )• -l-oo. Therefore, for 
any positive constants R and e one can find an integer I > 1 such that (j2.6p 
holds with Ci = • • • = C/ = and u — u{Q). Since K, contains the zero element, 
we see that Condition (H2) is satisfied. 

Example 2.3. Suppose that h is represented in the form h{t,x) — ho{Xt,x), 
where ho is a continuous 1-periodic function of time with range in V such that 
Jq ho{s) ds — 0. We claim that, for sufficiently large integers A > 0, Condi- 
tion (H2) is satisfied, provided that IC contains the zero element. 

Indeed, let us represent a solution of Eq. (|l.ip with f = h in the form 
u = w + g, where g{t) = h{s) ds. Then the function w must satisfy the 
equations 

vb + {w + g,V){w + g) - i^Aw + Vp = - A.g, div w = 0. (2.9) 
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The condition imposed on h imply that 5 is a l^-valued 1-periodic function such 
that 

sup ||g(i) ||y — > as A — > 00. 

Combining this with an argument similar to that used in Example 12. 21 one can 
prove that Eq. (|2.9I) has a unique 1-periodic solution w, which is globally expo- 
nentially stable as t -> 00. It follows that u = w + g is a globally exponentially 
stable solution for problem (ll.ip . (|1.2I) with f = h. As in Example 12.21 we 
conclude that the function u — u{0) satisfies the required property. 

2.2 General criterion for mixing and application 

In this section, we outline the proof of Theorem l2.11 which is based on two key 
ingredi ents: a coupling approach developed in |KSOH IMat02( IKPS02[ IMY02| 
IHai02[ IShi04[ in the context of stochastic PDE's and a property of stabilisation 
to a non-stationary solution of Navier-Stokes equations |BRS11| . We first recall 
an abstract result established in |Shi08| . 

Let X be a compact metric space with a metric dx and let (ufc, P„), u G X,he 
a family of Markov chains in X. We denote by Pk{u,r) its transition function 
and by and the corresponding Markov semigroups. Let (ni.,P„) be 
another family of Markov chains in the extended phase space X = X x X such 
that 

U,Pk{u,-)^Pk{u,-), KPk{u,-)^Pk{u',-) ioT u^{u,u') e X, k>0, 

(2.10) 

where Pk{u, F) denotes the transition function for {uk, Pu) and 11, 11' : X — > X 
stand for the natural projections to the components of a vector u ~ {u,u'). 
In other words, relations (|2.10p mean that, for any integer fc > 1, the random 
variable considered under the law P^, with u = (u, u') is a coupling for the 
pair of measures {Pk{u, ■), Pk{u' , ■)). We shall say that (iij,,P„) satisfies the 
mixing hypothesis if there is a closed subset B G X and positive constants C 
and (3 such that the following properties hold. 

Recurrence. Let t{B) be the first hitting time of the set B: 

t{B) = min{fc > : itfe e B}. 

Then t{B) is P^-almost surely finite for any u ^ X , and there are positive 
constants Ci and such that 

E„exp((5ix(B)) < Ci for tt 6 X. (2.11) 

Exponential squeezing. Let us set 

cr = min{fc > : dx{uk,u'k) > C e''^''}. (2.12) 
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Then there are positive constants C2, S2, and ^3 such that, for any u E B, 
we have 



P„{a = oo}>53, (2.13) 
E„(/{.<oo} exp((52a)) < C2. (2.14) 

The foUowing proposition is a particular case of a more general resuh estab- 
lishecHin |Shi08] (see Theorem 2.3). 

Proposition 2.4. Let (ufc,Pu) be a family of Markov chains for which there 
exists another Markov family (ttfc,P„) in the extended space X that satisfies 
relation (|2.10p and the mixing hypothesis. Then {uk,^u) has a unique stationary 
distribution fi G V{X), and there are positive constants C and 7 such that 

\\Pk{u,-)- ^i\\l<Ce-''^ forueX, k>0. (2.15) 

To prove Theorem I2.1[ we first observe that the Markov family (ufc,P„) 
possesses a compact absorbing invariant set X C H, and it suffices to study 
its restriction to X, for which we retain the same notation. We shah prove 
that (ufc,P„) satisfies the hypotheses of ProDOsition l2.4l A crucial point of our 
construction is the following result, which says, roughly speaking, that if two 
points u,u' G X are sufficiently close, then the pair (Pi(u, •), Pi(u', •)) admits 
a coupling whose components are close with high probability; cf. Lemma 3.3 
in [K50T] . 



Proposition 2.5. Under the hypotheses of Theorem \2.1[ there exists a constant 
d > such that for any points u,u' X satisfying the inequality \\u — u'\\ < d 
the pair (Pi(w, •), Pi(u', •)) admits a coupling {V{u,u'),V'{u,u')) such that 

P{\\V{u,u')-V'{u,u')\\ > l\\u-u'\\} <C\\u~u'\\, (2.16) 

where C > is a constant not depending on u,u' £ X . 

The proof of this proposition is based on a controllability property for the 
Navier-Stokes system and application of a concept of optimal coupling; see 
Sections 151 and [5Tl We now define a coupling operator 72. = {JZ,TZ') by the 
relation 

y [S[u,C),S[u' X )) for ||u - u'll > d, 

where C, and C,' are independent random variables whose law coincides with that 
of 771. Without loss of generality, we can assume that C, and C,' are defined on 
the same probability space as V and V . The required Markov family (Mfc,P„) 
is constructed by iterations of 7?.. Namely, let (flfc, Pfc), fc > 1, be countably 



^In IShiOSI . the proof is carried out in the particular case when B = B X B; however, the 
same argument applies in the general situation. 
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many copies of the probability space on which TZ is defined and let (H, J^, P) be 
the direct product of these spaces. We set 

Uo^{u,u'), Uk ='JZ{uk-i.,ujk), k>l. (2-18) 

The recurrence property will follows from approximate controllability (see Hy- 
pothesis (H2) in Section l2.ip . while the exponential squeezing will be implied 
by Proposition [531 

3 Control problem 

3.1 Squeezing 

In this section, we consider the controlled Navier-Stokes system (jl.ip - (|1.3p on 
the time interval Ji = [0, 1]. We assume that the function h belongs to the space 
H^{Di) (where Di — Ji x D) and denote by {ipj} the sequence of functions 
entering Hypothesis (HI) of Section 12.11 Extending the functions ipj by zero 
outside Q, we may regard them as elements of Hq{Di). We denote by 
the vector span of "01, . . . ,iprn endowed with the norm and by Bfj the ball 
in H^{Di) of radius R centred at origin. 

Theorem 3.1. Under the above hypotheses, for any R > and q £ (0, 1) there 
is an integer m > I, positive constants d and C, and a continuous mapping 

<P:BrxBh{R)^C{H,£,^), {h,uo) ^ v, 

such that the following properties hold. 

Contraction. For any functions h G Bfi and uo,uq G Bh{R) satisfying the 
inequality 

\\uo-uo\\<d, (3.1) 

we have 

\\S{uq, h) - S(uo,h + <P{h,uo){uo - uo))|| < 9 ||uo - ""oil- (3.2) 

Regularity. The mapping <P is infinitely smooth in the Frechet sense. 

Lipschitz continuity. The mapping <1> is Lipschitz continuous with the con- 
stant C. That is, 

\\$ihi,ui) - ^ih2,U2)\\^ < C {\\hi ~ h^Wn^ + \\ui - U2\\), (3.3) 

where \\ ■ \\c stands for the norm in the space C{H,£„i). 

Taking this result for granted, let us prove Theorem A stated in the In- 
troduction. Let us fix positive constants p and a and take an initial function 
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Uo € Bh{p). The boundedness of the resolving operator for the Navier-Stokes 
system imphes that the corresponding solution u satisfies the inequality 

\\u{t)\\<R foraUt>0, (3.4) 

where i? is a constant depending only on p. Let g G (0, 1) be such that e~" — q. 
Denote by d, C, and m the parameters constructed in Theorem 13. II Let us take 
any uq G H satisfying the inequality ||uo — uo|| < d and construct a control 
function 77 consecutively on the intervals Jk, k > 1. To this end, we denote 
by hk the restriction of the function h{t — k + 1) to Di and set 

r]{t) = ^(uq, hi){uQ - iio) for < t < 1. 

Then, by Theorem 13. 11 we have 

||m(1) - m(1)|| < q\\uo -uo\\ = e""||uo - uo||, 

where u{t) denotes the solution of p.ip - p.3p . (ll.Sp on the interval Ji. Assume 
we have constructed a control 77 on the interval (0, /c) with A: > 1, and the 
corresponding solution u satisfies the inequality 

\\u{l)'u{l)\\ \\uo-uo\\ (3.5) 

for Z = 1, . . . , fc. In this case, setting 

ri{k + t)^^{u{k),hk+i){u{k)-u{k)) forO<t<l, (3.6) 

we sec that (|3.5|) remains valid for ^ = fc + 1. Thus, wc can construct 77 on the 
half- line M-|_, and the corresponding solution will satisfy (|3.5p for all integers 
/ > 1. Combining this with the Lipschitz continuity of the resolving operator 
of the Navier-Stokes system, we see that ()1.6|) holds. Inequality (|1.7|) is a 
straightforward consequence of p.6p and (|3.5p . Finally, it is not difficult to 
check that if 77^ and 77^ are two controls corresponding to (uq,Uq,/i*), i = 1,2, 
then 

ll^fe - vIWlHd,) ^ ^1 (ll«o - "oil + ll"o - ^oll + max \\hl - /i?||ffi(D,)) > 

where Ci > depends only on p and a, rf^ stands for the restriction of 77* to Di, 
and h\ are defined in a similar way. This completes the proof of Theorem A. 

We now turn to the proof of Theorem 13.11 As was mentioned in the Intro- 
duction, a weaker version of this result was established in the paper jBRSll] . 
and its proof repeats essentially the argument used there. However, since the 
finite-dimensionality in time for the control and the regularity and Lipschitz 
properties of are important for the stochastic part of this work, we present a 
rather complete proof of Theorem l3.1l We begin with a description of the main 
steps and give the details in the next two subsections. 
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Step 1: Reduction to the linearised problem. Denote by u(t,x) the solution 
of (jl.ip . (|1.3I) issued from uq and corresponding to 77 = 0. In view of the 
regularising property of the Navier-Stokes system, for any interval J = 1) 
with S > we have 

ii e L^{J,H^ nV), dtiie L'^{J,V), (3.7) 

and the corresponding norms are bounded by a constant depending only on S 
and R. In particular, the truncated observability inequality holds for the adjoint 
of the Navier-Stokes system linearised around u; see Section [521 

A solution with a non-zero control rj is sought in the form u = ii + v. Then v 
must be a solution of the problem 

w + (w, V)t; + (u, V)t; + (w, V)u - i^Aw + = 7;(i,a;), divu = 0, (3.8) 
v\g^=0, v{0)^vo, (3.9) 

where vq = uq — uq- Together with Eq. (I3.8p . consider its linearisation around 
zero: 

V + {u,\/)v + {v,\/)u - i^Av + \7p = r]{t,x), divt; = 0. (3.10) 

Suppose that we have constructed r/ e £m such that the solution w{t, x) of p.lOp . 
(|3.9p satisfies the inequalities 

imi)ll<fl^oll, MU.<CM. (3.11) 

A standard perturbative argument shows that if ||do|| is sufficiently small, then 
the solution of p.8p . p.9p satisfies the inequality ||v(l)|| < 9||i^o||, whence it fol- 
lows that p.2p holds. Thus, it suffices to construct a continuous linear operator 
<?(/i, iio) : H £„i such that the solution w £ ^"1 of problem (|3.10p . (|3.9p with 
77 = <P(h,UQ)vQ satisfies inequalities p. 111) . 

Step 2: Application of the Foia§-Prodi property. Let {cj} be an orthonor- 
mal basis in H formed of the eigenfunctions of the Stokes operator L ~ —HA, 
where 11 stands for the Leray projection in L'^{D,'U.'^) (onto the closed sub- 
space H), let {ttj} be the corresponding (non-decreasing) sequence of eigenval- 
ues for L, and let Hn be the orthogonal projection in H on the vector space 
spanned by ei, . . . ,eN- Denote by TZ" : H x L'^{Di) Xi a linear operator 
that takes {vq , rf) to the solution w of p.lOp , p.9p and by 7?." its restriction to 
the time t. 

Suppose that for any integer > 1 and any 5 > we have constructed an 
integer m > 1 and a family of linear operators ^ = ^(^, wo) '■ H — > which is 
a Lipschitz function of its arguments and is such that 

\\UNnUvo,^ih,UQ)vo)\\ < C2S\\vol Mh,uo)\\c < C2, (3.12) 
where C2 > is a constant not depending on A^ and S. In this case, the Poincare 
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inequality and the regularising property of p. 101) imply that 

Wnfivo, <P{h, uo)vo)\\ = \\{I- IlN)nUvo,$ih, uo)vo)\\ + C26 \\vo\\ 

< a^Yf||7^^(^;o,^'(/i,wo)«o)||l + ^2^ Ikoll 

< C3a^Yi'(^o|| + \mh,uo)vo\\L^iD,))+C25\\vo\\ 

< {C3{C2 + l)a]^'lf + C2S)\\vo\\. 

Choosing N sufficiently large and S sufficiently small, we obtain the first in- 
equality in p. lip . The second is an immediate consequence of the continuity 
of TZ" and the boundedness of ^. 

Step 3: Minimisation problem. The construction of <P is based on a study of 
a minimisation problem for solutions of (I3.10p with a cost functional penalising 
the term ||nArii;(l)||. Namely, let us consider the following problem. 

Problem 3.2. Given a constant J > 0, an integer iV > 1, and functions vq G H 
and u E Xi satisfying (13. 7|) . minimise the functional 

J{w,C)^lj\mfdt+^\\U^wil)f 

over the set of functions (w,C) S -^i x L^(Z?i,M^) satisfying the equations 
w+{u,V)w + {w,\7)u-vAw + \7p = xiP,nO, divw = 0, u;(0) = uq, (3.13) 

where p = p{t, x) is a distribution in Di and P„j stands for the orthogonal 
projection in L^{Di) onto 

We shall show that Problem 13.21 has a unique solution {w, (), which satisfies 
the inequality 

i||nwi)f + IICIIi^p,)^^3lkof, (3.14) 

where C3 > is a constant not depending on N and 6. This will imply the 
required inequalities p.l2p . in which <Pvo = x(PmC)- Further analysis shows 
that the mapping u 1-^ ( is smooth from Xi to L'^{Di) and uniformly Lipschitz 
continuous on bounded balls. Since u is an analytic function of (/i, uo), this will 
complete the proof of the theorem. 

3.2 Minimisation problem 

Step 1: Existence, uniqueness, and linearity. Let us prove that Problem l3.2l has 
a unique optimal solution {w,() in the space X :— Xi x L'^{Di). Indeed, the 
function J : X — ;> R is non-negative and therefore has an infimum on any affine 
subspace of X. Denote by X^^ the affine subspace of X defined by p. 131) and 
by J* the infimum of J on Xyg. Let (u;„, Cn) S Xy^ be an arbitrary minimising 
sequence. Then {Cn} is a bounded sequence in i^(_Di), and without loss of 
generality we can assume that it converges weakly to a limit C. It follows that 
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the sequence of solutions Wn G Xi of problem p.l3p with C = Cn converges in 
the space A'l to a limit which satisfies p.lSp . The lower semi-continuity of 
the norm of a Hilbert space now implies that 

J{w, C) < liminf J(w„, («) = J*- 

Recalling the definition of J*, we conclude that J{w, C) = J*. 

To prove the uniqueness, note that any affine subspace is a convex set in X. 
Combining this property with the strict convexity of the norm of a Hilbert space, 
we see that if {wi,Ci), * = 1,2, are two optimal solutions, then Ci = C2. Since 
the solution of problem (I3.13P with a given ( e L^{Di) is unique, we conclude 
that wi = W2- 

Finally, it is a standard fact of the optimisation theory that the unique 
minimum of a quadratic functional under a linear constraint can be expressed as 
a linear function of the problem data. In the case under study, the corresponding 
operator depends on the reference solution u. We shall denote by !f'(it) the linear 
operator that takes uo to where {w, () is the optimal solution for Problem l3.2l 

Step 2: Regularity and Lipschitz continuity. We now prove that the map- 
ping regarded as an application from H to C{H, L^{Di)) is infinitely dif- 
ferentiable and uniformly Lipschitz continuous on balls. Let us denote by 
Wo = Wq{u,vq) G Xi the solution of problem p.l3p with ^ = 0. The linear 
constraint of Problem 13.21 is equivalent to the relation 

w = wo+7e"(0,x(PrnC))- (3.15) 

Setting A{u)( — 7?."(0, x(PmC)), see that (w,C) G ^ is a solution of Prob- 
lem [3?2] if and only if w is the global minimum of the function F : L^{Di) — >■ M. 
defined by 

F{C; u, vo) - i WCm'dt + 1 \\Ur,{wo + A(u)C) f . 

Using standard methods of the theory of 2D Navier-Stokes equations, we can 
prove that A{u) is an analytic function from Xi to C{H) which is uniformly 
Lipschitz continuous on bounded sets; e.g., see [Kuk82 or Chapter 1 in |VF88j . 
It follows that F satisfies the hypotheses of ProDOsition [?3] in which U = L^{Di) 
and Y = Xi X H. Thus, the unique minimum ^ of F is a smooth function 
of (m, uo) valued in L^{Di), and it is Lipschitz continuous on bounded subsets. 
Recalling relation p.lSp . we conclude that the unique solution (w, C) G X of 
Problem 13.21 is a smooth function on Xi x H which is Lipschitz continuous on 
bounded subsets. Since (w,C) linearly depends on vq, it is straightforward to 
derive the required properties of ^. 

Step 3: A priori estimate. From now on, we fix i?, (5 > and assume that 
ii G By^ J [R). We claim that there is an integer m > 1 depending on R and N, 
and a constant C > depending only on R such that, if w G By^ -^{R), then 
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inequality p.l4p holds. Indeed, note that the constraint given by p.l3p is 
equivalent to the equations 

w + iyLw + B{u,w)+B{w,u)=n{xiPmO), w{0) = vq, (3.16) 

where we set B{ui,U2) = n((iti, V)m2)- Thus, the pair (w,C) G X constructed 
in Step 1 is the minimiser of J under constraint (I3.14p . Applying the Kuhn- 
Tucker theorem (see Chapter I in |IT79| ). we can find functions \ E H and 
9 6 L^(Ji, F) such that, for any (r, G X, we have 

^\C, Odt + -^ {liNw{l), r(l)) + (A, r(0)) 

+ ( {e,f + i'Lr + B{u,r) + B{r,u)~xiPmO)dt^0. 
Jo 

It follows that 

e ~ vLO - B*{u)e ^0, (3.17) 
6(1) - -|nww;(l), C = PrnixO), (3.18) 

where B*{u) denotes the (formal) adjoint of the operator B{u, ■) + _B(-, ii) in i?, 
and Eq. p.l7p holds in the sense of distributions. Multiplying Eq. p.l7p by w 
and the first equation in p.l6p by 6*, adding together the resulting relations, 
and integrating over Di, we derive 

l\\IiMw{l)f + j'^ \m\?dt = -(0(0), z;o), (3.19) 

where we used relations p.lSp and the initial condition in p.l6p . Now note 
that p. 171) is equivalent to the backward Navier-Stokes equations (|5.12p . There- 
fore, by Proposition [521 the function 9 must satisfy the truncated observability 
inequality (I5.13p . Combining it with p.l9p . we obtain p.l4p . We have thus 
shown that (cf. ((3?T2t ) 

||^w7^^(^;o,xPm(^'(^)«o))|| <C5\\v4, ||3'(u)«o||l^(do <C\\v4. (3.20) 
3.3 Proof of Theorem EH] 

Let us define (l>{h,U()) as the linear operator taking vq to x(Pmtf'(u)wo), where 
u ^ Xi is the solution of problem (II. ip . (jl.6p with f — h and uq — uq. Standard 
regularity results for 2D Navier-Stokes equations imply that ii £ yg.i for any 
6 G (0, 1) (see Theorem 3.5 in Chapter III of |Tem790 . so that ^{ii) is weU 
defined. Moreover, the norm llullj;^! is bounded by a constant depending only 
on R and d. We claim that <?(/i, uq) satisfies the contraction property stated 
in the theorem. If this assertion is proved, then the regularity and Lipschitz 
continuity of 'F combined with similar properties of the resolving operator for 
the 2D Navier-Stokes system will imply the remaining assertions on <P. 
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Inequalities p.20p imply that the solution {w, () of Probleni l3.2l satisfies p.l4p . 
Therefore, choosing N and sufficiently large, we ensure that the func- 
tion w = TZ^{vo,(P{h,uo)vo) satisfies p. lip . Let us represent a solution of the 
non-hnear problem p.8p . (I3.9P with the right-hand side 77 = <l>(h,ua)[uQ — -uq) 
in the form v = w + z. Then z G Xi must be a solution of the problem 

i + (z, V)z + {it + w, V)z + (z, V)(u -f w) - v/\z + Vp = -{w, V)w, (3.21) 
divz = 0, 2|go=0, z(0) = 0. (3.22) 

Taking the scalar product of (|3.2ip with 2z in H and carrying out some standard 
transformations, we derive 

9t||2||2 + 2i.||Vzf < Ci(||zi + i«||i||z|| + ||z/;||i||w||)||Vz||. 

It follows that 

atlkf + HlVzlP < C^dlzi + u;||?|lzf + |lz.||?|hf ). 

Applying the Gronwall inequality and using the initial condition in p.22p . we 
obtain 




Recalling that is bounded by a constant depending only on R and using 

the second inequality in (|3.1ip . we derive 

sup ||z(l)|| < C^{R)\\w\W < Ci{R)\\uo - 7lof . 

0<t<l 

Since ||uo — wo|| < choosing 0? > sufhciently small ensures that 

lk(l)||<|lko-iiol|. 

Combining this with the first inequality of p. lip in which vq = uq ~ uq, we 
get (|3.2I) . This completes the proof of the theorem. 

4 Proof of Theorem 12.11 

We first outline the main steps. A well-known dissipativity argument shows 
that the random dynamical system defined by (j2.ip has a compact invariant 
absorbing set AT C H. Therefore it suffices to prove the uniqueness of an 
invariant measure and the property of exponential mixing for the restriction 
of (ufc, P„) to X. This will be done with the help of Proposition l2.4l Namely, we 
shall prove Proposition l2.5l and use relations (|2.16p and (|2.17p to define a Markov 
chain (tifc,P„) in the extended phase space X = X x X . This Markov chain is 
an extension of {uk, Pu) and possesses the recurrence and exponential squeezing 
properties of Section 12.21 and therefore the hypotheses of Proposition 12.41 are 
satisfied. This will complete the proof of Theorem 12. II 
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4.1 Reduction to a compact phase space 

Well-known properties of the resolving operator for the Navier-Stokes system 
imply that S satisfies the inequality 

\\Siuj)\\ <k\\u\\ + CiWfWmD,) for u e BhIR), f e l^{d,) 

where m: < 1 and Ci > are some universal constants. Let r > be so large 
ll'* + ^fc||L2(£)j) < r almost surely. Then, with probability 1, we have 

\\Siu,h + 7^k)\\ <^P + Cir hvu&BHip). (4.1) 

It follows that if i? > fr^j then the ball Bh{R) is invariant for the Markov 
chain {uk,Pu)- Let us take R = and denote by X the image of the set 

Bh{R) X ^L2(Di)(^) under the mapping S. Then X is an invariant subset for 
(u/c,Pm), and the regularising property of the Navier-Stokes dynamics implies 
that X is compact in H . Iterating (|4.1|) . we see that 

"^uiuk e X for fc > A:o(p)} = 1 for any u £ Bh{p), 

where ko{p) = (Inp + C2)/lnx-~^ with a large constant C2 > 0. It follows 
that (MfcjPu) has at least one stationary measure p, and any such measure is 
supported by X. It is easy to see that to prove (|2.8p . it suffices to estabhsh 
inequality (|2.15p . The latter is proved in the next three subsections. 

4.2 Proof of Proposition 12.51 

We shall apply Proposition 15.31 to construct a measurable coupling {V,V') and 
PropositionlOto prove ^A^. Namely, fix i? > so large that X C Bh{R-1) 
and \\h + ?7fc l|_f/i(Di) < i? — 1 almost surely. Denote by C, d > and m > 1 the 
parameters constructed in Theorem 13.11 with q — 1/4, define the Polish space 

Z = {{u, u) eX X X -.Wu- u'W < d}, 

and introduce the function e{u, u') ~ ^11"" ~ ^^'11 on the space Z . Let us consider 
the pair of measures (Pi(u, •), ■))■ By Proposition 15.31 with 6* = i, there 

is a probability space (51,7^, P) and measurable functions V,V' : il x Z ^ H 
such that for any {u,u') 6 Z the pair {V{u,u'),V' {u,u')) is a coupling for 
and 

P{\\V{u,u')~V'{u,u')\\ > <q|„-„-||/4(PiK-),AK,-))- (4.2) 

In view of (j5.6l) , to estimate the right-hand side of this inequality, it suffices to 
bound the function -ftr||n_„'||/4(-Pi(w, ■),Pi{u', ■)). 

To this end, we shall apply Propositions 15. 2[ 15.61 and Theorem 13.11 Let us 
endow the ball Bji C H^{Di) with the law A of the random variables 7/^,. Then 
Pi{u, ■) is the image of A under the mapping 77 1-^ S{u, h + rj) acting from Bj^ 
to H. Let £jn be the subspace entering Theorem 13.11 and let 

^■.BRxBHiR)^C{H,£m) 
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be a smooth C-Lipschitz mapping such that 

\\S{u,h + Tj) - S{u',h + r] + <P{h + T],u){u-u'))\\ < ^\\u-u'\\ (4.3) 

for any 77 S B;^ and any u,u' £ Bh{R) satisfying the inequahty ||w — < d. 
The existence of such mapping was established in Theorem 13. II Let us define a 
transformation >F = 'Fu,u' of the space H^{Di) by the relation 

<f (?y) = ?7 + x{^h + Ti\\i)^{h + ry, u){u - u'), 

where % : M+ — > is a smooth function such that x(r) = 1 for r < i? — 1 and 
x(r) = for r > i?. The choice of the constant R and inequality (|4.3p imply 
that 

||S'(u,/i + r/) - S'(w',/i + *„,„'(^))|| < for A-a.e. ?? e i?i(i:'i). 

Therefore, by Proposition [521 '^e have 

^||«-n'||/4(i^l(",-),^l(^',-)) < 2||A-<f',(A)||var. (4.4) 

Now note that the mapping ^ satisfies the hypotheses of Proposition 15.61 with 
a constant x proportional to ||m — Combining (|4.2p . ()4.4p . and ()5.22p . we 
arrive at the required inequality (|2.16p . 



4.3 Recurrence 

Let us recall that the probability space (r2,jF, P) and the X-valued Markov 
chain {tt^ = (m^jM^)} on it were defined in Section We shaU denote by P„ 
the probability measure associated with the initial condition u. Let us set 

B ^{u^ (u, u') e X : \\u - It'll < d}, 

where d > is a small constant to be chosen later. In this and next subsections, 
we shall prove that the recurrence and exponential squeezing properties are 
satisfied for (ttfc,P„) with the above choice of B. 

To prove that t{B) is P„-almost surely finite and satisfies p.lip . it suffices 
to show that 

supP,{T(B) < 1, (4.5) 

where £ > 1 is an integer. Indeed, if this inequality is proved, then using the 
Markov property, for any integer to > 1 we can write 

Fu{r{B) > ml} = E„P„{r(S) > ml \ 

= E„(/,(b)>(„_i),P„(„^_^,Jt(S) > I}) 
<pPu{TiB) > (m-l)£}. 



19 



where {J-k} stands for the filtration generated by the Markov family (Mfc,P„). 
By iteration, the above inequality implies that 

P„{r(B) > me} < for all m > 1. (4.6) 

A simple application of the Borel-Cantelli lemma now implies t{B) is P„-almost 
surely finite. Furthermore, inequality (j4.6l) immediately implies that t{B) sat- 
isfies (EHI) . 

We now prove (I4.5p . Let 77 and r/ be the random variables entering the 
definition of the coupling operator IZ (see (|2.17p ) and let 

Cfc(w) = ri{ujk), C'ki^) = V'i^k), k>l. 

Then {Cfc, C^i fc > 1} is a family of i.i.d. random variables in L^{Di) whose law 
coincides with that of rjk- Let u ^ H he the point defined in Hypothesis (H2). 
Since X is a closed absorbing subset, we have u € X. Using the definitions of TZ 
and X, we can write 

Pu{t{B) > £} = Pu{t{B) > e, \\ue ~ u',\\ > d} 

= Fu{r{B) > i, \\Sfiu, Ci, . . . , Q) ~ Se{u, ([,..., CM > d} 

< n\\se{u, Ci, . . . , Ci) - SfXu\ Ci, . . . , c;)ii > 

= l-F{\\SfXu,Ci,--.,Ci)~Siiu'X'i,---Xi)\\<d} 

<i-p{e,df, 

where the operator Si is defined in Hypothesis (H2), and we set 
P(£, d) - inf ¥{\\S,{v, Cu...Xd-u\\< d/2}. 

Thus, inequality (|4.5p will be established if we show that P{e, d) > for a 
sufficiently large integer £. 

To this end, let us denote by i the integer defined in Hypothesis (H2) with 
e = d/A and i? > so large that X C Bh{R)- We fix u g X and denote by 
Cj. = C^ G /C some vectors for which (|2.6p holds. Let (5 > be so small that for 
any C,k G /C, fc = 1, . . . satisfying the inequalities — || < we have 

||5,(«,Ci,...,6)-"ll <rf/2. 

Define the event T{v,t) = {||Cfe - Cfcll < (5, A: = 1, . • . , ^} C fJ. What has been 
said implies that 

\\St{vXi,---Xd~u\\<d/2 foru; Gr(t;,^). 

We see that 

P{i,d) > inf P(r(z;,^)) > inf J] HWCk ~ QW < S} 

> (mf P{||Ci-^l| <<5})', (4.7) 
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where we used the fact that Cfc are i.i.d. random variables. It remains to note 
that the function ^ i— IP{||Ci ~ Cll < ^} is lower semicontinuous and strictly 
positive on A^, and therefore the right-hand side of (|4.7|) is a positive constant. 
This completes the verification of the recurrence property. 

4.4 Exponential squeezing 

Let us take arbitrary initial points G X such that {u,u') G B. Then, 

by (Pl^ and (PTT)) . we have 

Pu{\\ui - u[\\ <^\\u~ u'W} >l-C\\u- u'\\. (4.8) 

For any integer n > 1, let us set 

r„ = {\\uk - 411 < hW^k-i - 4-1 II for 1 < /c < n}. 

Note that ||ufc - u'J < 2-'=||mo - MqH for 1 < fc < n on the set r„. Combining 
this observation with inequality ()4.8|) and the Markov property, we derive 



Pn(r„) = E„(/r„_,P„{||u„ - <|| < i||M„_i - 4_i|| I Tn-i}) 
-E„(/r„_,P„„_,{||iii -411 < i||Mo-4ll}) 
> (l-2i-"C||u-u'||)P„(r„_i). 

Iterating this inequality, for any n > 1 we obtain 

n 

P«(r„) > l[{l-2^-''C\\u~u'\\) >l~C4u~u'\\. (4.9) 

k=l 

Let us define the Markov time a by relation (|2.12l) with C = d and /3 = In 2. 
Then, for any u B, with Pu-probability 1 we have 

{cr = cx)} = {||ufe - 411 < 2-''d for A: > 1} 

CXD 

D {\\uk - 411 < ihk-i - 4-ill for fe > 1} = fl r„, 

n=l 

whence it follows that 

P„{a = cx)}= lim P„(r„) > 1-Ci||u-m'||. 

Choosing d > sufficiently small, we arrive at (|2.13p . 
To prove inequality (|2.14l) . it suffices to show that 

Pu{<T = n} < C22"" for n>l, ue B, (4.10) 

where the positive constant C does not depend on m. To this end, note that 

{a = n} = {a > n-l}n {\\u„ - u'J > 2-''d}. 
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Using the Markov property, inequality (|4.8p . and the fact that — mJ.|| < 2 ''d 
on the set {a > k}, we derive 

PuW - n} = ^u{l{a>n-l}Pu{\\Un ~ u'J > 2""^! 

= IE«(/{.>„-i}P«„_i{||tii - u'lll > 2-"d}) 

< IEu(/{a>„-l}Pu,._i{||ui - u[\\ > i||uo - u'o\\}) 

< 2^-''CdFu{<J >n-l}. 

We thus obtain inequaUty ()4.10p with C2 = Cd. The proof of Theorem 12.11 is 
complete. 



5 Appendix 

5.1 Optimal coupling 

Let X be a Polish space with a metric d and let ^i,/i2 G 'P(^). Recall that 
a pair of X-valued random variables (^1,^2) is called a coupling for {^1,^2) if 
= — 1, 2. We denote by n(^i, ^2) the set of all couplings for {^1,^2)- 
Let us fix £ > and define a symmetric function d^ : X x X ^ M.hy the relation 

Jl if d(Mi,W2) > e, 
I II d(ui,U2) < e. 

Definition 5.1. We shall say that a coupling (^1,^2) € n(/zi,/X2) is e-optimal if 
it minimizes the function (Ci, ^2) ^ Ec?e(Ci, C2) defined on the set of all couplings 
for {iJLi,iJL2)- That is, 

IE4(a,6) = , inf E4(Ci,C2). (5.2) 

In particular, for e = we obtain the usual concept of maximal coupling of 
measures; e.g., see |ThoOO) . 

Let us denote by Ce(/ii,/X2) the quantity defined by the right-hand side 
of (j5.2l) and call it the e-optimal cost for the pair (/xi,/i2). Note that, in view 
of (|5.ip and (|5.2I) . for any e-optimal coupling (^1,^2) we have 

PM6,6)>e} = a(Mi,A^2). (5.3) 

Thus, to estimate the probability of the event that the distance between the 
components of an e-optimal coupling for {pi,^2) is larger than e it suffices to 
estimate the corresponding e-optimal cost. We now establish a simple result 
that enables one to do it. Let us introduce the following function on the space 
ViX) X P(X): 

Ke{fJ.l,ti2) = sup((/,/ii) - {g,ti2)), (5.4) 
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where the supremum is taken over ah functions /, g G Cb{X) satisfying the 
inequahty 

/(wi) - 5(1*2) < 4(^1, M2) forui,U2eX. (5.5) 
Then, by the Kantorovich duahty (see Theorem 5.10 in [VilOQj ). we have 

Ke{fJ.l,f^2) ^ Ce{fll,fl2) (5.6) 

Thus, to estimate the e-optimal cost it sufhces to estimate the function K^. The 
foUowing proposition reduces this question to a "control" problem. 

Proposition 5.2. Let X be a compact metric space with metric d, let Ui,U2 
be two X -valued random variables defined on a probability space {fl,J-,F), and 
let /ii,/i2 be their laws. Suppose there is a measurable mapping if" : 17 — > such 
that 

d{Ui{uj), [/2(!f'(w))) < £ for almost every uj e fl, (5.7) 
where e > is a constant. Then 

ifs(/^l,M2) <2||P-<f,(P)||var. (5.8) 

In particular, any e-optimal coupling (^1,^2) for the pair (/ii,/i2) satisfies the 
inequality 

FM6,6) >£} <2||P-<Z'*(P)||var. (5.9) 

Proof. Inequahty (|5.9p is a straightforward consequence of (|5.8p , (15.61) , and (|5.3p . 
To prove (|5.8p . we use an argument apphed in the proof of Lemma 11.8.6 
in |Dud02] . Namely, note that if /, g G Cb{X) are such that (|5.5p holds, then 
the function h{u) = sup„g^(/(w) — de{u,v)) satisfies the inequalities h < g and 

/(mi) - /i(u2) < de(Mi,M2), \h{ui) - h{u2)\ < I forui,U2GX. 

It follows that 

(/,Mi)-(5,/^2)=IE(/(C/i)-g(t/2)) 

< E (/(C/i) - h{U2 o<Ir))+E {h{U2 o <F) - h{U2)) 

< E4([/i, U2 0^) + {hoU2, M^)) -{hoU2,P) 

< \{hoU2,F) - {hoU2,MP))\. 

It remains to take the supremum over /, g and to note that right-hand side of 
this inequality does not exceed that of (|5.8p . □ 

We now study the question of existence and measurability of an e-optimal 
coupling. Let {Z, B) be a measurable space and let {/if , 2 G Z}, i = 1,2, be two 
families of probability measures on X such that the mapping z 1-^ fif is measur- 
able from Z to the space 'P(^) endowed with the topology of weak convergence 
and the corresponding tr-algebra. The following proposition establishes the ex- 
istence of an "almost" e-optimal coupling that is a measurable function of z, 
provided that X is a subset of a Banach space. 
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Proposition 5.3. In addition to the above hypotheses, assume that X is a 
compact subset of a Banach space Y . Then for any positive measurable func- 
tion e{z), z ^ Z, and any 9 G (0,1) there is a probability space (r2,J^, P) and 
measurable functions (w) : il x Z — > y, i = 1, 2, such that, for any z ^ Z , the 
pair (^f , ^1 ) is a coupling for (fif , /i| ) and 

^d,^,^ie^,e2)<Ceei.M,^i'2)■ (5.10) 

Proof. The existence of a measurable optimal e(z)-coupling would be an imme- 
diate consequence of Corollary 5.22 in |Vil09| if the function d^(z){ui,U2) was 
continuous and independent of z. Since this is not the case, we now outline the 
modifications that are needed to prove the existence of a measurable coupling 
satisfying (|5.10p . 

Let us write Z as the union of countably many disjoint measurable subset Zk 
such that the image of the restriction of e{z) to any of them is contained in a 
bounded interval separated from zero. For instance, we can take Zk — {z G Z : 
{k + 1)^^ < £{z) < k^^} with fc > 1 and use a similar partition for (1, +oo). If 
we construct measurable couplings on each set Z^ for which (|5.10p holds with 
Z = Zk, then by gluing them together, we get the required random variables , 
i = 1,2. 

To construct a measurable coupling on Zk, we first reduce the problem to 
the case when e = 1. Let us consider stretched measures /if defined by the 
formula /if(r) = ^^{e{z)T) for F G By- The measures /if are supported by a 
compact set Xk C Y, and if (CiiCI) is a measurable coupling such that 

Kdi{eiA2) < CeiK.f^l) for any z e Zk, 

then £(2)(^f,^|) is a measurable coupling for the original measures that satis- 
fies (|5.10p . Thus, we can assume from the very beginning that e = 1. 

Let d : X X X — ^ K be an arbitrary continuous symmetric function such that 

di{ui,U2) < d{ui,U2) < dg{ui,U2) for all ui,U2 6 X. (5.11) 

By Corollary 5.22 in |Vil09| . there is a probability space (51, T, P) and measur- 
able functions S^f : Z x ^ X such that (ff,f|) is a coupling for (/^f,/i|) for 
any z ^ Z , and 

^d{£,l,Q)^ ini Ed(Ci,C2). 

(Ci,C2)en(^if 

Combining this relation with (|5.11l) . we arrive at (|5.10p with e = 1. □ 



5.2 Truncated observability inequality 

Let 5 £ (0, 1) and let ii G ys^i be an arbitrary function. Consider the adjoint 
equation for the Navier-Stokes system linearised around u: 

g+ {u,V)g+ {Vg,u) + i^Ag + VTT = 0, div.g = 0. (5.12) 

Let us fix an open set Q C {6, 1) x D, an orthonormal basis {^j} C L^(Q,M^), 
and a function x G Co°(Q) and denote by Pm the orthogonal projection in 
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L^((5,R^) onto the m-dimensional subspace spanned by (pj, j = 1, . . . ,771. Re- 
call that stands for the vector span of the first N eigenfunctions of the Stokes 
operator. The following result is a simple consequence of the observability in- 
equality for the linearised Navier-Stokes system (see |FE991 UmaO 1[ lFTJGIP04j ) . 

Proposition 5.4. For any S,p > and any integer N > 1 there is an integer 
m > 1 such that if u G By^ ^ (p) , then any solution g ^ Xi of Eq. (|5.12p with 
g{l) G Hn satisfies the inequality 

||5(0)|| <C||P™(x.9)|L.(^^), (5.13) 

where C > is a constant depending only on 6 and p. 

Proof. We essentially repeat the argument used in the paper |BRS11| in which 
inequality (15.131) is proved for a time-independent function x ^-^d & projec- 
tion P„j acting in the space variables. We claim that if g S is a solution 
of ([51^ with g{l) g H^, then 

where Ci > is a constant depending on S, p, and N . Once this inequality is 
established, the required result can be derived by using the fact that 

||(/-Pm)f|| <Sm\\v\\i for any v e iJi(Q,K2), 

where {dm} is a sequence going to zero as m — ?> 00; cf. proof of Proposition 5.3 
in IBRSllj . 

Suppose that (|5.14p is false. Then there are functions m„ G By^ ^ (p) and 
solutions gn G Xi of ()5.12|) with u — Un such that 

5„(l)eiJA,, ||5„(1)|| = 1, (5.15) 
IU5n||^i(^^) >n||xff„||i.(^^). (5.16) 

Now note that ||(7„(1)||^2 < C2 for all n > 1, whence it follows, by standard 
estimates for the 2D Navier-Stokes system, that 

hnWxt + \\9n\\L'^(JM-i) + || ^tffn || L2 ( j y) < C3, U > 1. 

Passing to a subsequence, we can assume that {{*„} and {gn} converge weakly 
(in appropriate functional spaces) to some functions ii and g, respectively, such 
that u e By^^{p), g G Afi is a solution of (|5.12|) . and ||5(1)|| = 1. Moreover, it 
follows from (|5.16p that xg = 0- Let us show that the latter is impossible. 
Indeed, let an interval (a, 6) C J and a ball B C D he such that 

X{t, x) > a> for (i, x) G (a, b) x B. 

Then, by the observability inequality (see Lemma 1 in jFCGIP0"4] ). we obtain 

||ff(a)|| < C'4||g||L2((a,h)x_D) < C4a-^\\xg\\L^Di) = 0. 

The backward uniqueness for the linearised Navier-Stokes system (see Sec- 
tion II.8 in [BV92]) implies that g{t) = for a < i < 1. This contradicts 
the fact that ||5(1)|1 = 1. The proof is complete. □ 
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5.3 Minima of second-order polynomials on Hilbert spaces 

Let U and Y be real separable Hilbert spaces and let F : J7 x F — > K. be a 
function of the form 

y) = {QyU, u)u + {ay, u)u + by, 

where Qy G C-{U) is a self-adjoint operator, ay, by e U for ah y and (•, ■)u 
stands for the scalar product in U. We assume that 

[Qyu, u)u > c\\u\\lr for all ueU, ye Y. (5.17) 

In this case, it is easy to see that for any y e Y the function u i-> F{u, y) 
has a unique global minimum u* = u*{y). The following simple proposition 
establishes some properties of u*. 

Proposition 5.5. Let us assume that the functions 

y Qy, y ^ {ay, by) (5.18) 

are infinitely smooth from Y to the spaces C{U) and U x U , respectively. Then 
the unique minimum u*{y) is an infinitely smooth function ofyGY. Moreover, 
this implication remains true if the property of infinite smoothness is replaced 
by Lipschitz continuity on bounded balls. 

Proof. Inequality (|5.17p implies that u*{y) is the only solution of the linear 
equation 

{duF){u,y)^Q. (5.19) 

Therefore, the required smoothness of u* will be established if we show that 
the implicit function theorem can be applied to (|5.19l) . This is a straightfor- 
ward consequence of the fact that duF is an infinitely smooth function of its 
arguments, and its derivative in u coincides with 2Qy. 

We now prove that if Oy, by, and Qy are Lipschitz continuous on bounded 
balls, then so is u*{y). Indeed, fix i? > and take two points yi,y2 G By{R). 
It follows from (|5.19l) and the explicit form of the derivative duF that 

Qvi{u*{yi) -u*{y2)) = {Qy, - Qy,)u*{y2) + ^{oy, - ay,). 

Taking the scalar product of this equation with u*{yi) — u*{y2) and using in- 
equality (I5.17P and the Lipschitz continuity Qy and Oy, we obtain the required 
result. □ 

5.4 Image of measures under finite-dimensional transfor- 
mations 

Let X be a separable Banach space endowed with a norm || • || and represented 
as the direct sum 

X = E + F, (5.20) 
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where E C X is a finite-dimensional subspace. Denote by P^; and Pp the pro- 
jections corresponding to decomposition (|5.20|) . Let A G "PiX) be a measure 
that can be written as the tensor product of its projections \e = (P_b)*A and 
A_F = (Pi?)*A. Assume also that A has a bounded support and that Xe possesses 
a C^-smooth density p with respect to the Lebesgue measure on E. The follow- 
ing simple result gives an estimate for the total variation distance between the 
measure A and its image under a diffeomorphism of X acting only along E. 

Proposition 5.6. Under the above hypotheses, let ^ : X ^ X he a transfor- 
mation that can he written in the form 

<I'{u) =u + 'P{u), ueX, 

where <P is a -smooth map such that the image of $ is contained in E and 

||'^(wi)|| < ||<?(ui) - '^(u2)|| < - U2II for all ui,U2 e X. (5.21) 

Then the total variation distance between A and its image under W admits the 
estimate 

\\X-M^)\Ur <Ck, (5.22) 
where C > is a constant not depending on k. 

This proposition is a simple particular case of more general results presented 
in Chapter 10 of |BoglO| . However, for the reader's convenience, we give a 
complete proof of Proposition l5.6l 

Proof. Inequality (|5.22|) needs to be proved for x ^1, because it is trivial if x 
is separated from zero. Furthermore, since A has a bounded support, there is 
no loss of generality in assuming that <P vanishes outside a large ball. 
We first note that if / € Cfc(X), then 



(/,A) = / Xpidw) / f{v + w)p{v)dv. 

J F J E 

Since v t-^ v + <P{v + w) is a C^-diffeomorphism of i? for |>ir| ^ 1, we see that 
(/,iZ/,(A))= [ Xpidw) [ f{v + w + <P{v + w))p{v)dv 



^ det (/ + (D^) ie^iv')+w)) 

where Q^{v') denotes the solution of the equation v -\- (^{v -\- w) = v' and 
stands for the Frechet derivative of ^. Combining the above formulas, we get 



S{f,X):={f,M^))-if,X)= / Xpidw) / fiv + w)Aiv,w)dv, (5.23) 

Jf J e 

where we set 

A/ N P(0u.(w)) , s 
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The >r-Lipschitz continuity of ^ implies that 

\det{l + {D$){Qyj{v) + w)y^ - l| = 0(^<) for all veE,weF, 

where 0(>f) denotes any function whose absolute value can be estimated by C>c 
(uniformly with respect to all other variables). It follows that 

det(/+ (e„(u) +w)) 
Furthermore, using inequalities (|5.2ip . we derive 

l|e.(z;)|| < i\\v\\ + \\QUv) ~ < {M + M). (5.25) 

Substituting (|5.24p into (|5.23p and using inequalities (|5.25p . we obtain 

|^(/,A)|< / XAdn.) I \fi^ + -^)\:rT7^f^fj^^ 

Jf Je det[I + {D<P){Q^{v) +w)) 

< II./II00 / Xpidw) [ \p{v)- p{v + <P{v + w))\dv + Oi>c). 
Jf Je 

Taking the supremum over / G Ch{X) satisfying ||/||oo < 1 and noting that the 
integrals on the right-hand side can be taken over sufficiently large balls, we 
obtain 

||A-tf'4A)|lvar < ^ J XF{dw) J \p{v)- p{v + ^{v + w))\dv + 0{}<) 

Bf{R) Bf{R) 

<^ [ XF{dw) [ [ \\/p{v + e<P{v + w))\dedv + 0{>c), 



2 

Bf{B.) Bf(,B.) 







where we used the first inequality in (|5.2ip . It is straightforward to see that the 
right-hand side of the above inequality can be estimated by Cx. □ 

Acknowledgment. The author is grateful to Lihu Xu for stimulating dis- 
cussion. 
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